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O; 

. In previous works, we have constructed a refined version of the Gribov-Zwanziger action in 4 dimen- 

sions, by taking into account a novel dynamical effect. In this paper, we explore the 3-dimensionaI 
case. Analogously as in 4 dimensions, we obtain a ghost propagator behaving like in the in- 
frared, while the gluon propagator reaches a finite nonvanishing value at zero momentum. Simulta- 
neously, a clear violation of positivity by the gluon propagator is also found. This behaviour of the 
propagators turns out be in agreement with the recent numerical simulations. 

I. INTRODUCTION 

m ' 

J> ' Lately, the infrared behavior of the gluon and the ghost propagator in SU{N) Yang-Mills theories has been 
exhaustively investigated by many research groups. The low energy behavior of these propagators are of great 
interest as they might provide helpful information on various aspects of color confinement, which is still far from 
being understood. Previous results on the gluon and ghost propagators in the Landau gauge have reported an 
enhanced behavior for the ghost and a suppressed gluon propagator vanishingat zero momentum. This behavior 
00 was supported by both numerical simulations (ll |2j] and analytical studies I|3t S 0, H, Ht llOl . 
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Nevertheless, recent lattice data on bigger volumes point towards a ghost pr opagato r which is no longer en- 
hanced and a gluon propagator which attains a finite value at zero momentum [12l, Il3l, [T3, [Tsll . Recently, several 
analytical approaches 1,16.,,17„,18.,.19„,2Q1 have been worked out, being in agreement with these new data. 



Among these approaches, we have settled for a framework within the Gribov-Zwanziger approach, successfully 
employed in the study of the propagators in 4D [17, 18] in the Landau gauge, which is defined as d^Api = 0. We 
recall here that this condition does not uniquely fix the gauge freedom, as there can be configurations A|,, gauge 
equivalent to A^, which also fulfill djiA'^ = fT]. The Gribov-Zwanziger action allows for a (partial) resolution of 
this problem of gauge (Gribov) copies in a local and renormalizable setting LZ'^^iSJ. 

Our refined framework was constructed as follows. We have added two extra terms to the ordinary Gribov- 
Zwanziger action, without destroying the renormalizability of the action. Let us already recall here that these addi- 
tional terms precisely correspond to extra, yet unexplored, dynamical effects associated with the Gribov-Zwanziger 
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action iITtIITsII . The first extra term corresponded to the introduction of a novel mass operator, while the second term 
represented an additional vacuum energy term which was required in order to remain within the Gribov region Q. 
This region is defined as the set of field configurations fulfilling the Landau gauge condition, 3^A^ = 0, and for 
which the Faddeev-Popov operator, 

M'^ = -d^ {df,S"^'+gr^A'f?) , (1) 

is strictly positive, namely 

Q = {a;, 8f,A« = 0, M'^ > 0} . (2) 

This region is bounded by the horizon, 8Q, where the first vanishing eigenvalue of M."^ appears. Doing so, a large 
number of gauge copies is already excluded, as their appearance is precisely related to the existence of zero modes 
for A^"* f?]. So far, we have only worked out the 4-dimensional case. The implementation of the Gribov-Zwanziger 
approach in 3D has not yet been carried out. The 3D case is also conceptually different from the 4D case due to the 
superrenormalizability, and there is no running of the coupling constant g-^. Hence, it is useful to study the 3D case 
in detail and make a comparison with the 4D case as well as with the available 3D lattice data. 

Also in lattice simulations, one has to deal with the existence of gauge copies in the Landau gauge when studying 
the propagators. The Landau gauge is numerically implemented by minimizing the functional J d^x over its 
gauge orbit. Notice that the Gribov region Q corresponds in fact to the set of all relative minima of J d^x A^. An 
ideal implementation of the Landau gauge would correspond to finding the absolute minima of that functional. 
Even numerically, this is an extremely hard task to achieve, let stand alone to do it in the continuum formulation. 
Nevertheless, the gauge field configurations employed in the numerical evaluation correspond to relative minima 
of the functional J d^x A^, so that they belong to the region Q. Let us also mention that, at present, it is unknown 
how to restrict in the continuum the path integration to gauge fields corresponding to absolute minima of / d^x A2. 
That set of absolute minima forms a subset of the Gribov region Ci, known as the fundamental modular region A, 
A C n. This means that Ci itself is still plagued by additional gauge copies |2ill22ll23ll . Albeit it has been argued 
in ll23l that these additional copies should not contribute to the correlation functions, so that there should be no 
difference between the regions D. and A, we point out that the restriction to O via the local and renormalizable 
Gribov-Zwanziger approach is the best one can do for now in the continuum. Summarizing, both in the lattice 
and continuum study of the gauge theory, one is looking at the same objects, e.g. gluon and ghost propagator (see 
Appendix C and Section V), in such a way that the issue of gauge copies has been at least partially handled. We 
refer to e.g. [12, 13, 14, 15, 25] for lattice works. 

For the benefit of the reader, we review the original construction of the Gribov-Zwanziger action in sec- 
tion II. Section III presents a detailed motivation of why we should include effects of a mass term like S^^ = 
Af^ J d^x {q)cp — cUco) to the Gribov-Zwanziger action Sqz- Also, the renormalizability of Sgz + %<p is discussed. 
The requirement of renormalizability/ consistency with the symmetries of the model puts rather severe restrictions 
on the possible additional operators which can be introduced in the theory. We also provide some details on the 
breaking of the BRST symmetry. In section IV we discuss the need of the inclusion of an extra vacuum term 

Sen = / d^x^^^^^gA^. Again, we prove the renormalizability of the total action Stot = Sgz + S-^<p + Sen- In 

section V, we investigate in detail the gluon and ghost propagator and already show that for 7^ 0, the ghost 
propagator is not enhanced and the gluon propagator is finite at zero momentum, in accordance with the lattice 
data of 1121 [isl. [Til. [Tsl] . In Section VI a variational technique is scrutinized in order to determine the value of M^. 
This mass parameter is thus not added by hand to the original Gribov-Zwanziger action, but is determined in a 
self -consistent way. Using this variational technique, we shall also investigate the one loop ghost propagator in the 
full momentum range. Finally, our conclusion is given in section VII. 

II. THE ORIGINAL GRIBOV-ZWANZIGER ACTION 

We summarize the construction of the Gribov-Zwanziger action in d dimensions, where it is understood that we 
actually work in dimensional regularization, with d = 3 — 2£. We start from the following action |j9|] 

Sstart = Sym + SLandau " lY J d'^xf^'A^ (M'^)"' A], , (3) 
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where 

Sym = I / d'^^f «,f;, (4) 

is the classical Yang-Mills action and 

SLandau = / d'*x(fc«a,,A« +c'd,Dfc'') (5) 

denotes the Landau gauge fixing and the ghost part. The part of expression (|3]l proportional to the Gribov mass 
parameter 7 is the nonlocal horizon function, which implements the restriction to the first Gribov region, with the 
proviso that this 7 is not free, but subject to the horizon condition |9|] 

{h{x)) =d{N^-l) , (6) 

where h{x) is the so called horizon function 

hix) = gY^'At (M-^X'f^'K ■ (7) 



This nonlocal horizon function has also received attention from the lattice community, see for example [26], where a 
lattice formulation of the horizon function can be found. 



In order to find a manageable local quantum field theory, one adds extra fields Vlir ^Iff localize the 

nonlocal part of the action (|3]l. Doing so, the Gribov-Zwanziger action becomes [[9l l27ll 

Sgz = So - / d'x {j'gf'Alcp'; + I'gf'A;-^'; + d[N'- l) 7^) , (8) 

with 

So = Sym + I d^'x (b''df,Al, + c''df,{D^,cy) 

' d'^X (f^dy {DyCp,Y - Widy {DyCV,)" - g {dyW^) /"^^ {DyC)'' ff) , (9) 



whereby (^l^^i (p'^^ ^re ^ pair of complex conjugate bosonic fields, whereas {^^f , '^"f^ ^re anticommuting ghost 

field. Based on a global U{f) symmetry, f = d (N^ — 1), with respect to the composite index i = {}i,c) of 

the additional fields {^'^^ , cp'^ ,Ld"f , Ll)'^^ we have introduced a more convenient notation {^"^ , cp"^ ,Ld"^ , cv'^f^ = 

{^1 , (pi itd^ , co") . Notice that for d = 3, dim(g^) = 1 and dim(7^) = 3/2. Defining the quantum effective action 
r by means of 

e"^ = / dYe-s , (10) 



where Y is a shorthand for all the fields, then it is easily shown that the horizon condition (|6j is in fact equivalent 
with the requirement that 
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2 - (11) 



is fulfilled for 7^ 7^ 0, or 



gr''A;{cp';+v';)) = -^d (n^ - 1 • (12) 



This means that 7 will become fixed in terms of the natural scale of the theory which, in 3D, is provided by the 
coupling itself. We thus expect to find 7^ g^- 

As it has been shown in 0,113, H^], the action ^ is renormalizable to all orders, and is thus suitable for quantum 
computations. 
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III. THE GRIBOV-ZWANZIGER ACTION COMPLEMENTED WITH A NEW OPERATOR 



A. Proposal 

As carried out in the 4-dimensional case fv?, ^T§], we introduce the local composite operator (LCO) into the 
action ||8ll. We shall couple this mass term to the action using a source /, and for renormalization purposes, we have 
to add this term in a BRST invariant way. Therefore, we shall see that the o^-sector must gain the same mass. Hence, 
we consider the following action: 

S^^ = J d'x[-]{Y^cp1-Wlco1)] , (13) 

with Sqz the original Gribov-Zwanziger action l(8|. The new source / has the dimension of a mass squared. There- 
fore, we shall often use the notation = /. 



B. Motivation 



Let us briefly repeat the motivation for adding an extra term to the Gribov-Zwanziger action. As pointed out in 
lIT^ , the fields {^"^, (p'^fi'^'^^, '^"i^ introduced to localize the horizon function appearing in ||3| are interacting fields 
corresponding in fact to the nonlocal dynamics associated with the horizon function. Therefore, these fields will 
develop their own quantum dynamics. We draw attention to the fact that the A- and cp- fields are intimately entan- 
gled, since there is a quadratic Aip-mixing term present in the tree level action (|8j, namely gy^ f^'^ A'^^{cp^^ + '^^^)- 
Hence, we expect that any effect in the ^-sector will immediately reflect in the gluon sector, altering, in particular, 
the behavior of the gluon propagator at zero momentum. This is nothing else than saying that the horizon strongly 
influences the gluon dynamics. 

On top of the previous argument, we can give a second argument why one should add the local composite operator 
(LCO) {^(p — coco) to the action. The horizon condition corresponds to 

gr''A;{cp';'+^'p) = -2d (14) 

i.e. to a condensate proportional to 7^. One might then also expect that {fcp — coco) will be nonvanishing too when 
7^ 7^ 0. Hence, one is almost obliged to incorporate the effects related to the operator '^cp — Tvco. To make this 
argument more explicit, let us compute the perturbative value of the condensate {fcp — coco) . We start from 



_ > dW{}) 

(W-^^>pert = 



(15) 



with W(/) the generating functional defined in our case as 

-W(j) ^ /[dY]e-S' (16) 



and with S' the extended Gribov-Zwanziger action given in flSl l. At lowest order, one finds 
Making use of 



(271)'' V q^ + J J J {2Tz)'i ' I ' ' < ' J (2n) 



^''^ . In(<j2 + col)+[^ He,' + cl) - / ln(^2 + (18) 



(27r)'* " ^' J {2n)^ " ~' J {In) 
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with 



(19) 



where = Ig^Ny'^, and after employing a standard integral in dimensional regtilarization. 



{Inf ' ' (Atz) 

we find the following ultraviolet /infte result 



^i'^(^' + '«') = -S;w2rM/2)' (20) 



W(/) = -3(n2 - 1)^ + f-a;3+ - a;^ + f/A . (21) 



2g2Af en 

Using this explicit expression, we can easily obtain the perturbative value of the condensate ((TSt , reading 

(^<p - ZUa;)pg,t = V2 A w 0.056(n2 - 1)A , (22) 

where A is the nonzero solution of = 0. Since at one loop 



2Ng2 2 ' '7 (271)" 



we find that 



with 



A=^^^N, (24) 

_ , N3(n2-1) 

10368714 (^^^ 

We notice that the one loop vacuum energy corresponding to the Gribov-Zwanziger action is positive. The same 
feature was also observed in the one loop 4D case |18]. 

Besides this perturbative value of the condensate (|22] l, it could also be possible that another, non-perturbative, 
value emerges. However, to calculate this non-perturbative value, one should calculate the Legendre transformation 
of the generating functional W(/) to obtain the effective action T(a), where the absolute minimum configuration 
^ {^(p — loco) would correspond to the true, ener geti cally favoured, vacuum. Here, we limit ourselves to 
observe that, similar to what was encoimtered in 4D, |18|] , no such value has been found at one loop, meaning 
that higher order contributions have to be taken into accoimt. For the interested reader, we have written down the 
corresponding details in Appendix lAl 



C. Renormalizability 

To prove the renormalizability, let us start from the following action, 

S" = S' + Slco / 
Slco = / d^W} , (26) 

with S' the action proposed in ( [131 , and with the LCO parameter p a new dimensionless quantity. This term in p 
is, in principle, needed to take into account potential divergences proportional to which are allowed by power 
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TABLE I: Quantum numbers of the fields. 



counting and by the symmetries of the action . This term would ensure multiplicative renormalizability of the 
functional However, although this term cannot be ruled out at the level of the algebraic analysis, an analogous 

proof as in 4D can be written down allowing us to consistently set p = fl^ . 

To study the renormalizability of the Gribov-Zwanziger action, it is highly useful to embed it in an extended 
action, which reduces to the original model in a specific limit |9]. Doing so, we may have a larger number of Ward 
identities at our disposal, which are powerful tools to construct the most general possible counterterm. Therefore, 
proceeding as in the 4D case 118.1 , we shall start with the following action 

^ = So + Ss + Sext + S-^ip + Slco / (27) 

with So given in ||8]l, S-^^p in ((TSl l, and 

Ss = sj d'x (-W; {D.cpiY - (D„aJ,)" - u^f) , 

= J d'^x (-M«/ {D.cpiY - gU;:/"'' (D.c)" cpj + Lf- (D.o;,)' 

Sext = / d'^x (^-K; (D.c)" + igL7'"'^c''c^) . (28) 

We introduced new sources M°/, V^' ,U", N", K", and L", which are necessary to analyze the renormalization of the 
corresponding composite field operators in a BRST invariant fashion. The BRST operator s is defined through 

sA; = - (D,,c)" ,sc" = ^gf^'c^c' ,sc" = b" ,sb" = ,s< = a;f ,sa;« = ,sa;« = ^« ,s^° = , 

sU^' = Mil , sM\l = , sVf = N«' , sN«' = , sX« = 0, sL" = , s/ = , (29) 

so that s is nilpotent, = 0. One can easily see that the action Z is indeed BRST invariant, sE = . We underline 
that the mass operator itself, — Zocv, is also BRST invariant. In TABLE I, we have summarized for all the fields and 
sources their mass dimension, ghost number and Q y^-charge, which is defined by means of the diagonal generator 
Ujj of the global JJ(/) symmetry. We have chosen these mass dimensions such that in any case, the action of the 
BRST transformation s raises the dimension by 1/2. 



At the end, we can give the sources the following physical values 



M 



W" 



phys 



phys 



yah 
^ liv 



phys 



^2 rafcv 



phys 



phys 



phys 



0, 



(30) 



in order to recover the physically relevant action S", given in | |26|| . If the action E is renormalizable for any value of 
the sources, it will of course be for the specific values l l30t . This is an example of the fact that it can be highly useful 
to use a slightly more general version than the original action, whereby a more powerful set of Ward identities can 
be invoked to prove e.g. the renormalizability of S'. We have collected the details of the algebraic renormalization 



^ We refer to >29i1 for more details concerning the LCO formahsm and the precise role of the LCO parameter 
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analysis in Appendix B. The main conclusion drawn is that the action E (|27t is multiplicatively renormalizable to all 
orders of perturbation theory. 

It is important to notice here that it is the Ward identities defining the extended Gribov-Zwanziger action (|27|l 
which dictate which terms can be added to the action, without jeopardizing the symmetry content in general and 
the renormalizability in particular. From this perspective, the mass term ((T3)l displays remarkable features. In fact, 
its introduction turns out to be compatible with both Ward identities and renormalizability, as originally proposed 
inllTJ. 

D. The breaking of the BRST symmetry 

We would like to draw attention to the fact that the actions S and S' are BRST invariant if 7^ = 0, since then 
the extra fields can be integrated out, and we are left with the original Yang-Mills action in the Landau gauge. 
Evidently, we then have {lpq> — aJcv) = — {s{aJ(p)) = 0. It is only when 7^ 7^ that we can have such a nonvanishing 
condensate, since in that case the BRST symmetry generated by s is no longer preserved, since sS' 7^ 0. More 
precisely, the "horizon terms" oc 7^ in equation (|8} are not BRST invariant. It is an important observation that 
the Gribov-Zwanziger action is not BRST invariant: the BRST symmetry is explicitly broken by soft terms ^ 7^^. 
Nevertheless, it is worth emphasizing that this breaking can be kept under control at the quantum level |18]. In 
particular, the introduction of the local sources U, V, M and N, see also [30], nicely allows one to embed the action 
into the larger BRST invariant action, thereby allowing one to prove the renormalizability. When the sources attain 
their physical values, eqs.JSOjl, the exact Slavnov-Taylor identities of the larger action induce the corresponding 
softly broken identities for the physical action S". The operators coupled to these sources are exactly those relevant 
for the discussion of the broken Slavnov-Taylor identity; see a similar discussion in JSOfl. 

The introduction of the horizon function thus explicitly breaks the BRST invariance of the Landau gauge fixed 
Yang-Mills action. In particular, in [18], it has been shown that the origin of this breaking is deeply related to the 
restriction to the Gribov region Q. More precisely, it turns out that infinitesimal gauge transformations of field 
configurations belonging to the Gribov region Q give rise to configurations lying outside of Q. The appearance of 
the BRST breaking looks thus rather natural, when we keep in mind that the BRST transformation of the gauge field 
is inherited from its infinitesimal gauge transformation. Let us elaborate a bit on this here. The starting point is that 
we must restrict the domain of path integration to the region Q. Let us consider the full set of fields present, Ap, c, c, 
b, and potential other fields. These fields can be seen as the coordinates of a space !F. We can define a manifold over 
J-', by means of the action functional, 

S : Ap,c,c,b,... -> S{Af„c,c,b, . . .) e R. (31) 

To be more precise, we must restrict ourselves to O, so we must consider the fimctional 

Srestricted : A^, C,C,b,...(z T\Ap G ^ S{A^„ C, C, &,...) G R . (32) 

In a pictorial way, we can imagine this as some kind of "cylinder" in with O as ground surface, and the other 
unrestricted fields (c,c, . . .) describing the height. Notice that we do not say what the action S precisely is, we only 
assume that it is BRST invariant. Now, we consider a particular point (A*, c* ,c* , b* , . . .) very close to the boundary 
of the cylinder, thus A* is located very close to the inner boundary of O, and 3^ A* = 0. We can decompose A* as 

a; = fl;, + C^, (33) 

with Cj, G dCl, thus C^, lies on the Gribov horizon. The shift fl^, is a small (infinitesimal) perturbation. Obviously, 
dpCp = d^a^t = 0. We then find 

= C,, +fl,,+D;,(C)a; + ... (34) 

a; 

for the gauge transformed field. Since G dCl and setting to equal to the zero mode corresponding to C^,, yields 



^^,A^, = 3,,D,,(C)a; = 



(35) 
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In addition, A^, also lies very close to the boundary 30, however on the other side, as shown by Gribov in f7]. Let us 
now consider the infinitesimal BRST shift of the coordinate set (A*, c* , . . .), which is given by 

a; ^ a; + 0D,,c* (36) 

for the gauge field, with 6 a Grassmann number. We did not specify yet the choice of our particular ghost coordinate. 
We take c* = with 6' another Grassmann number'^, and oj the zero mode. Since the BRST was supposed 

to be a symmetry of the cylinder, the transformed coordinate set of (A*,c*, . . .) should still be located within the 
cylinder. However, by construction of Ap and c* , we do end up outside of the cylinder. This contradiction means 
that maintaining the BRST symmetry is not possible when restricting to O. 

One can also imagine a configuration A** G O not located close the boundary 30, thus 9fiA** = and 
-d^i D,, ( A** ) > 0. If we are to assume that the BRST transformation of A**, 

A** + 6D^iC** , (37) 

with c** arbitrary, would remain within O, we are forced to conclude that d^D^[A**] [dc**) = 0, thus dc** would 
then be a zero mode, again in contradiction with the hypothesis that A** is not located on (or close to) the boundary 
30, meaning that there are no such zero modes. 

An interesting property worth being mentioned is that, despite the loss of the BRST symmetry, one can still 
use the associated broken Slavnov-Taylor identity to derive relations amongst several Green functions. This has 
exhaustively been studied in [18], and can be easily transported to the 3D case. We therefore refer the reader to 
Gi for any detail concerning the BRST breaking and its consequences. In particular, we have also algebraically 
motivated that the (controlled) BRST breaking allows the Gribov parameter 7 to become a physical parameter. 
If the BRST symmetry would be preserved, 7 would merely play a role akin to that of an unphysical gauge parameter. 

We emphasize here that the soft breaking of the BRST is introduced in such a way that one keeps the nilpotency of 
the BRST operator. This is an important point, as it enables one to make use of the notion of the BRST cohomology. 
Indeed, as discussed earlier in this section, the proof of the renormalizability of the action is exactly possible by 
making use of the BRST cohomology in the extended model | |27|| , which enjoys the BRST symmetry. In fact, one can 
define the physical operators in the Gribov-Zwanziger model as those obtained from the physical operators Ophys 
in the extended model l|27|l , upon taking the physical limit l|30|l of these operators Ophys/ which are nothing else than 
the cohomology classes of the nilpotent BRST operator. The latter are precisely given by the gauge invariant singlet 
operators constructed with the field strength and its covariant derivative. In other words, due to the form of the BRST 
operator, the physical operator content of the theory is left imchanged, being identified with the colorless gauge 
invariant operators. As a consequence, both gluons and ghosts are excluded from the physical spectrum. Moreover, 
due to the presence of the Gribov parameter 7 and of the mass M, the behavior of the correlation function of gauge 
invariant operators will also get modified in the infrared region, as it can be inferred from the expression of the 
resulting gluon propagator. For instance, although the gluon propagator turns out to exhibit positivity violations (see 
section V, eq. (|50|l ), one might expect that the correlation functions of gauge invariant operators, like e.g. (F^(x)F'^(i/) ), 
could display a real pole in momentum space, which would be related to the mass of a glueball bound state. This 
topic is currently under study. At asymptotically large momenta, one can neglect the soft BRST breaking term, in 
which case we are reduced to the normal Yang-Mills physical modes. The degrees of freedom corresponding to the 
fields cp, Ip, CO and cJ will decouple from the physical spectrum, as these fields form BRST doublets, see ||29] |, and it is 
well known that these become trivial in the BRST cohomology [31,] . 



E. A few more words about the intricacies of 3D gauge theories: infrared problems and ultraviolet finiteness 

In the previous subsection, we mentioned the case 7^ = 0. Strictly speaking, the 3D theory will not be well 
defined in that case. In the absence of an infrared regulator, the perturbation theory of a super-renormalizable 
3D gauge theory is ill-defined due to severe infrared instabilities 13211 . This can be intuitively understood as the 



^ Notice that 69' is a normal number, thus we can divide by it. 



+ _ _ 



FIG. 1: The one loop corrected ghost propagator. 

coupling constant carries the dimension of a mass. In the absence of an infrared regulator, the effective expansion 
parameter will look like /p with p a certain (combination of) external momentum /-a. For p ^ g^, very good 
ultraviolet behaviour is apparent, but for p <^ g^, infrared problems emerge. The presence of (a) dynamical mass 
scale(s) ni <x g^ could ensure a sensible perturbation series, even for small p, as a natural expansion parameter is 
then provided by g^/ m. From this perspective, a nonvanishing Gribov mass 7^ could also serve as infrared cut-off. 
This feature is also explicitly seen from eq. l|24|l , from which an effective dimensionless expansion parameter can be 
derived as 

- ' -0.6, (38) 



(47r)3/2A 2V2^ 

a quantity which is at least smaller than 1. The inverse factor (in)^^^ is the generic loop integration factor generated 
in 3D. 



Although this is not the main concern of this paper, it might be interesting to perform higher order computations 
to effectively find out whether all infrared divergences are absent when 7^ 7^ 0. 

In fact, it is possible to couple the regulating mass term ^rrP-A^ to the Gribov-Zwanziger action. We did not con- 
sider this in the current paper, but in 4D this has been discussed in full detail in [18]. All Ward identities and relations 
between renormalization constants are maintained. Moreover, the form of the propagators is only quantitatively in- 
fluenced by this additional mass rn^, whereby the main consequences of the mass related to '^(p — Tvco are preserved 
(see later) [18]. In the presence of ^rn^A^, all Z-f actors are 1, or said otherwise, there are no ultraviolet divergences 
when computing Green functions I33l l33l . Having a look at the relations | |B18|| and l lB19t , this also means that any 
other Z-factor is 1, and hence the Gribov-Zwanziger theory is completely ultraviolet finite, including the vacuum 
functional, since there is no independent renormalization for it: the potential divergences related to 7^ are killed 
by the already available Z-factors, which are themselves trivial, and we already know that there are no divergences 
related to g^]. It is understood that, if needed, the mass related to is brought back to zero, as the other mass 
parameters are expected to cure the theory in the infrared. This should be checked case by case. For the purposes of 
this paper, we shall later see that everything works out fine without a mass coupled to A^. 



IV. THE GRIBOV-ZWANZIGER ACTION SUPPLEMENTED BY AN EXTRA VACUUM TERM 



A. Proposal 

We propose to add an extra vacuum term to the action | [T3t , i.e 



where the pref actor 2^^^^== is chosen for later convenience. 



d'x g 72/ . (39) 



B. Motivation 



We shall now explain the need for the inclusion of this vacuum term. So far, we have altered the original Gribov- 
Zwanziger action by adding a mass operator to the action. However, we have to be careful that with this addition 
we are not leaving the Gribov region when performing calculations. Initially, staying inside the Gribov horizon was 
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assured by the horizon condition ((TT] |. This condition is equivalent to demanding that the Faddeev-Popov operator 
defined in |[l) is positive, i.e. — 3^,D"/' > 0. In turn, looking at the Landau gauge fixing ||5]l, this is equivalent to 
demanding the positivity of the inverse ghost propagator, which can be deduced from FIG.[TJ giving. 



with 



J«*'l(l+(7(fc2))+0(g^ 





1 








N2 - 1 y 



N 



N2 - 1 fc2 



A^q {k-q)^ky I 
(27r)4 {k-qf " 



(27r)4 {k-qf \ 



(40) 



(41) 



and with ^A",Af,y the gluon propagator (see later). We can now rewrite the ghost propagator by performing a 
resummation as 

1 1 



k^ 1 - Cr(fc2) 



O(g^) 



(42) 



as we are only working up to order g^. This corresponds to the usual resummation of a set of connected diagrams 
into the inverse of the one loop IPI ghost self energy. We can therefore write the inverse ghost propagator as. 



(43) 



g-\k') = k'{l-a{k'))+0{g''), 
Demanding the positivity of the previous expression is translated in the so called Gribov no-pole condition 

(T{k^) < 1 , (44) 

which is the key-point of this investigation. Without the inclusion of this extra vacuum term, we shall demonstrate 
that is impossible to satisfy the no-pole condition. We recall that this no-pole condition was the basis of the original 
Gribov paper [7]. 

We have introduced a new parameter g which is still free and needs to be determined. We notice that this vacuum 
term is, in a sense, comparable to the vacuum term already present in the action ||8]l, i.e. — / d'^xii (n2 — l) 7^. In fact, 
in the original Gribov-Zwanziger formulation, this term was also necessary to stay within the horizon. Analogously 
as 7 is fixed by a gap equation, we shall introduce a second gap equation to determine Q. Proceeding as in the 4D 
case, we shall impose that 



d(T{0) 



M2=0 



0, 



(45) 



which assures a smooth limit to the original Gribov-Zwanziger case when 

C. Renormalizability 



0. 



The renormalizability of the following action, 

Stot = S + Sen 

= Sgz + / d'^x 



(^«<-aJ«a;«)+2^i^£=ilg7' 



(46) 



can be proven analogously as in 1181 . The vacuum term Sen will not give rise to any additional counterterms, hence, 

eo7o/o _ G7^/ 



and consequently, adding this extra term does not give rise to a new renormalization factor. 



(47) 



(48) 
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V. THE GLUON AND GHOST PROPAGATOR 

We shall use the following conventions for the gluon and the ghost propagator, 

{a;{-p)aI{p)) = 3"'vip^)r,M , {c^i-p)c\p)) = 3"'g{p') , m 

where Vfivip) = &}iv — is the transverse projector. 



A. The gluon propagator 

The tree level gluon propagator corresponding to the action | [26l l is given by 

with 

= Ig^N'y^ . (51) 
This particular propagator l|50] l displays the following properties already at tree level: 

• ©(p^) is infrared suppressed due to the presence of the mass scales and A*. 

• T>{0) = i.e. the gluon propagator does not vanish at zero momentum if is different from zero. 

These properties seem to be in qualitative accordance with the lattice data llT2l[T3l[T3l . We also want to stress that 
the mass term related to '^cp — uJw plays a crucial role in having T>{0) 7^ 0, since in the standard Gribov-Zwanziger 
scenario, the gluon propagator necessarily goes to zero. 

It is instructive to determine the one loop value of the gluon propagator near zero momentum in order, for 
example, to produce a numerical estimate to compare with other methods. Again this calculation is similar to that 
performed in 4D and we record those features which are different for our 3-dimensional analysis. First, we note 
that for this exercise we will follow [35] where the form of the gluon and (p field 2-point mixing term is defined 
differently. For this intermediate calculation here it is appropriate to use those conventions and method of llssll 
since they have been demonstrated to be consistent with the gluon suppression and ghost enhancement of the 
original Gribov analysis at two loops in MS in 4D. At the end we have been careful in converting back to the 
main conventions of this article in deriving the one loop freezing value of the gluon propagator. Our method 
involves computing the matrix of 2-point functions comprising the gluon- field sector at one loop in the MS 
scheme and then inverting this at one loop to determine the quantum corrections to the tree propagators. However, 
as we are ultimately only interested in the zero momentum limit we restrict ourselves to evaluating the 2-point 
functions in the zero momentum limit from the outset using the vacuum bubble expansion. In this all fourteen 
contributing Feynman diagrams are expanded in powers of the external momentum but truncated at 0((p^)^). 
Unlike in 4D all the diagrams are ultraviolet finite and the renormalization of all the parameters is trivial. In other 
words one simply replaces bare quantities by their corresponding renormalized ones. In effect this is a trivial MS 
renormalization. In addition the basic one loop vacuum bubble integral is simple to compute and is given by 
differentiating l l20ll with respect to m^. To determine the freezing value of the gluon propagator we have adapted the 
computer programme used for the 4-dimensional computation, fl^, written in the symbolic manipulation language 
Form, [36], to the 3-dimensional case. As indicated already this is the main reason for adopting the conventions 
of [35] here and it requires the replacement of the basic one loop 4-dimensional integrals by their 3-dimensional 
counterparts in the FORM programme. Whilst this may seem a trivial exercise there is an important aspect to note. 
In performing automatic Fe5mman diagram computations, where the graphs are generated electronically with the 
QGRAF package, [37], one has always to ensure the correctness of the result. Paradoxically an ultraviolet divergent 
4-dimensional calculation is easier to check than a finite 3-dimensional one. The reason for this resides in the fact 
that for the former the correct non-trivial renormalization constants have to emerge, and in the case of a gauge 
theory, these have to satisfy the Slavnov-Taylor identities. In a finite calculation the luxury of this check is absent. 
However, in the present case the adaptation of a verified programme with minor changes gives us confidence in 
the eventual value we will derive. Moreover, a reasonable degree of consistency with other results, such as lattice 
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methods, adds to our confidence in the result as will become evident later. 

For completeness, the Landau gauge propagators we use are, for an arbitrary colour group, 

(AUv)A''(-v)) - ^""(P^ + M^) 



where f^'^ are the colour group structure constants and the appearance of 1 / \/2 derives from the conventions of 
iSSl . We formally define the matrix of one loop corrections to the 2-point functions as 

y_ry2jrcab _(-p2_^^2)^flc^brf J^[^^j:cab Qgac^bd ^ jrace jrbde ^ j^jabe jrcde ^ globed ^ J ' W^/ 

in the basis where the common tensor Pj,^(p) has been removed, d^*"^'* is defined by, Issll . 

and (r^)bc: = — zy*"^ is the adjoint representation of the color group generators. Given the set of formal one loop 
2-point functions, it is easy to invert the matrix to one loop and formally derive the corresponding matrix of propa- 
gators, as 

(p^+M^) rep rcpq 

7^^ fpcd 1 xcpKdq 1 7^ fcdr fpar 

[{p^f+M^p'+CAi^V (p^+M^) ^ Y+M^)[Yf+MY+CA'y^V ^ 

^jrpcd Q^cp ^dq _^ J jrcpe jrdqe _^ ^jrcde jrpqe _^ J S ~^ ) ■ (55) 

We refer to Section 3 in [35] for a detailed discussion about the color structures emerging in this matrix. 

As we are specifically interested in the gluon propagator the only quantity of importance here is A and it is 
formally the same as in fl^. In other words. 



A 



[(p2)2 + M2p2 + C^74]2 



Nf yx - CajHn + !i) (p^ + + CaI* (Q + CaR + W) .(56) 



This concludes the derivation of the formal aspects of the computation of the one loop propagator corrections. The 
actual values of the 2-point function contributions now need to be inserted from the vacuum bubble expansion. 

Let us point out here that we are looking at the connected gluon two-point function, which is the relevant quantity, 
also measured on the lattice. The reader might notice that the lowest order part of the gluon propagator in l l55t can 
also be obtained by integrating out the additional [ip, cp)- fields in the quadratic approximation and construct the 
tree level gluon propagator in this fashion. In the case = 0, this also corresponds with the result for the gluon 
propagator from the original semiclassical approach by Gribov It is worth mentioning that this particular kind 
of propagator has been used frequently as an ansatz for a long time to do lattice fits, see e.g. |l2j]. 

A key difference in evaluating the relevant fourteen Feynman diagrams derives from the dimensionality of the 
basic d-dimensional integral defined by 

for a generic mass fi. On dimensional grounds has mass dimensions of {d — 2). In four dimensions this of course 
corresponds to a mass dimension of two. By contrast in three dimensions this drops to one. However, in each of the 
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tree propagators there is a denominator factor which is quadratic in which gives two roots. In 4D this leads to 
a very simple partial fraction into two terms and hence the simple evaluation of the basic vacuum integral In 
3D the situation is more complicated. The same partial fraction emerges but one has to choose the sign of the square 
root of the mass term appearing as fi^ for the evaluation of the integral in terms of an object of mass dimension one. 
In principle this could lead to several analyses depending on the choice of sign in the square root. Another way of 
viewing this is to realize that in three dimensions one has to consider the common propagator factor as quartic in 
y^p^ and find four roots. For completeness these are 



m 



m 



4 



In the limit where M ^ one recovers the four more recognizable poles of the usual Gribov propagator 



lim mi 

M-»0 



lim mi 



(58) 



(59) 



To resolve which is the correct choice of signs for the integral, we recall that the gap equation for the Gribov mass 
derives from a one loop calculation which also involves the basic integral Examining that calculation in order 
to have a non- trivial solution with a positive coupling constant, one has to take and mi. Interestingly in our 
calculation the only square root which remains in the leading order vacuum bubble expansion relevant for the gluon 
propagator freezing is \/ IVP- + I^/CaI^ which is always real. So there are no issues concerning the relative sizes of 
and which could have led to a complex value in the square roots of | [58l l. Given these considerations and being 
careful to revert to our general conventions, we finally find. 



17, 



VM2 + 2A2 13,4VM2 + 2A2 



^A" 



M4 - 4A4 



-A^M^ 
3 



+ 2A2 10 8 VM2 + 2A2 7 
(M4 - 4A4)2 + 3 (M4 - 4A4)2 + 4 



-2A2 



(60) 



From the previous expression, we can still conclude that for a non zero value of M^, "0(0) 7^ 0. 



B. The ghost propagator 

We have found that the one loop corrected ghost propagator can be written as 

1 1 



fc2 1 - (7(fc2) ' 

where cr(A:2) is the following momentum dependent function 

'^(^ ) = ^ ilnf (ic-,)2,4 + M2,2 + A4 ^^-('^) • 

Calculating this integral explicitly, we find 

Ng2 r 1 



Zlk^n I m2 - VM4 - 4A4 V VM4 - 4A4 

-k'^n+^ (ik^ + lvf- Vm4 - 4A4^ ^ arctan 



V^A:Vm2 - Vm^ - 4A4 ^(m^- Vm^ - 4A4 

^/2k 



(61) 



(62) 



3/2 



Vm2 - VM4 - 4A4 



1 



M2 + VM4 - 4A4 V VM4 - 4A4 



M2 



^ V2A:3 \/m2 + Vm4 - 4A4 - (m^ + Vm4 - 4A4 



-fc^TT + i (2^2 + + Vm4 - 4A4) ^ arctan 



Vm2 + VM4 - 4A4 



3/2 



(63) 
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In order to find the behavior of the ghost propagator near zero momentum we take the limit fc^ ^ in equation 

-(0) = M /^l + = ^' + ^' (64) 

which can of course be obtained by taking the hmit A:^ — > in expression | [651 l. Similarly, one can check that — > 
for fc2 ^ 00 and/ or oo. 

Before drawing any conclusions, we still need to have a look at the gap equations, which shall fix as a function 

of M2. 



C. The gap equations 

We begin with the first gap equation fTTjl in order to express A as a fimction of M^. The effective action at one loop 
order is given by 

T? = -diN^ - 1)7^ + I'-i^, ,^M- + - 1) / A 1- + + . (65) 

^ ^ ' ^/Ig^ 2 ^ 'J {2nf q^ + M^ 

With A* = 2g^Ny^, we rewrite the previous expression. 



and apply the gap equation l(TT)l , 



In 3 dimensions, the integral in this gap equation is finite, resulting in, 

67r VM2 + 2A2 

This expression will fix A^ as a function of M^, i.e. A2(m2) once we have found an explicit value for g. 
This explicit value for g will be provided by the second gap equation l(45)l . From expression | |62|| , one finds 

cr(0) = g^N- ' ' ' 



d J {2nY + MY + \^ 



.d-1 f d'q 1 , ,,2„2..d-l f d^^ 1 



d J {InYq^ + M^q^ + X^^^^^ d J {2tzY q^ q^ + + \^ ' ^^^^ 

Therefore, we can rewrite the first gap equation | |67|| as 

= .(0)-1-„^1^,^n/^^^,-J^+4^ ,70, 

The second gap equation can then subsequently be obtained by acting with on the previous expression and 
setting = 0. Doing so, we find 

d-1 [ d'^q 1 1 1 
d ^ J {2nyqW + X'^{0)^'^A^0) " ^ ^ 



by keeping (|45ll in mind. Setting = in |[68] | yields. 



2A2(0) = ^ . (72) 

OTI 
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Proceeding with equation (|7TJ, we find the following simple solution for g, 

1 2 v'^N 
127T V2 A(0) 

In summary, the following expression, 

1 ^ (74) 

fixes A2(m2). 

D. The ghost propagator at zero momentum 

At this point, we have all the information we need to take a closer look at the ghost propagator at zero momentum. 
From equation l|64t and l(74)l . we find 

^'°)^(^-)(-^)--^- 

From this expression we can make several observations. Firstly, when = 0, we find that (7(0) = 1, which is 
exactly the result obtained in the original Gribov-Zwanziger action Consequently, the ghost propagator (|6l]l 

is enhanced and behaves like 1/A:'^ in the low momentum region. Secondly, for any > 0, u(0) is smaller than 
1. By contrast, without the inclusion of the extra vacuum term, (7{0) would always be bigger than 1, which can be 
observed from expression ((TOl l. Therefore, it is absolutely necessary to include this term. With (t{0) smaller than 1, 
the ghost propagator is not enhanced and behaves asl/k^. 

VI. STUDY OF THE DYNAMICAL EFFECTS RELATED TO 

A. Variational perturbation theory 

In this section, we shall rely on variational perturbation theory in order to find a dynamical value for the hitherto 
arbitrary mass parameter Afi. Specifically, we follow fisllsgll and introduce as a variational parameter into the 
theory by replacing the action (|46l l by 



{l-£'')M^j 



(^«<-a;?a;«)+2^i^5ilg7' 



(76) 



where £ serves as the loop counting parameter, and formally £ = 1 at the end [18]. In this fashion, it is clear that 
the original starting action Sqz has not been changed. We have in fact added the terms in M^, and subtracted them 
again at k orders higher in the loop expansion. We shall set A: = 2 as we are working up to one loop. Taking k = 1 
would destroy the effect of the vacuum term Sen/ which would be inconsistent as explained before. 

We shall discuss a possible option to fix M^. For this, we can expand any quantity Q^"\ which is evaluated to 
a certain order n, in powers of £, cut this series to the order n and set £ = 1. We can then require that = 0, 
i.e. the principle of minimal sensitivity jioll . This latter requirement can be well motivated since an exact calculation 
of Q using the action l(76] |, taking the full nonperturbative nature of the theory into account, would lead to an M^- 
independent result, since £ = 1 after all. At a finite order in £ however, some residual dependence will enter the 
result. In this way, we can hope to capture some relevant nontrivial information yet at finite order, encoded in the 
parameter M^. We then mimic the independence on of the exact result just by imposing = and thereby 

fixing M^. If, however, no optimal value for is found, one can impose that ^g^^f^ = li4Q1 . 
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FIG. 2: 17(0) in function of in units ^ = 1 



fc2 



0.05 0.1 0.15 0.2 0.25 0.30 0.35 0.40 0.45 0.50 0.55 



0.19 0.16 0.14 0.12 0.10 0.08 0.06 0.04 0.03 0.02 0.01 



TABLE II: Some Mi- for different fc^ in units €^ = 1. 



B. The ghost propagator 



In order to obtain a dynamical value for the ghost propagator in the presence of M^, we shall now utilize the 
variational procedure outlined in the previous subsection. We start with the expression ( |62l l for u(fc^) and replace 
g-^ with £g-^ and with (1 — ^^)M^, expand up to order £ and set £ = 1. Doing so, we recover again the same 
expression ([62J. Following an analogous procedure for the gap equations, also results in the same expression |[74|l . 
This latter equation determines A^(M^), which can be plugged in the expression of cr[k^), thereby making cr[k'^) 
only a function of M^, next to the momentum dependence. 

Firstly, let us investigate cr{k^) at zero momentum, which is the key-point of this paper. FIG. |2] displays cr[Q) as a 
function of M^. We observe that i7(0) is indeed smaller than 1 for all > as already shown analytically in the 
previous section. We also find a smooth limit of c(0) for — > required by the second gap equation ||45] |, as can 
be seen from the left figure. According to the principle of minimal sensitivity, we have to search for an extremum. 



I.e. 



3t7(0) 
3M2 



3V(0) 



0. Unfortunately, there is no such an extremum present. Nevertheless, we do find a point of inflection, 

2 



^ = at 



0.185 ( ^ 



. Taking this value for M , we find 



(7(0) 



0.94, 



(77) 



which is indeed smaller than one and results in a non-enhanced behavior of the ghost propa gato r. This value needs 
to be compared with the lattice value of (7(0) = 0.79, which can be extracted from the data in Ual . 

Secondly, let us have a look at this point of inflection, when "turning on" the momentum k^. As shown in TABLE 
HH we observe that will decrease, until it will vanish at k^ w 0.55. The corresponding cr(fc^) is displayed in FIG.|3l 
We thus see that we find some kind of a momentum dependent effective mass M^{k^), which disappears when k 
grows. This could have been anticipated, as we naturally expect that the deep ultraviolet sector should be hardly 
affected. Let us also notice here that cr[k^) is a decreasing function, as can be explicitly checked from FIG.|3] This of 
course means that we are staying within the horizon for any value of the momentum. 

To compare our results with available lattice data, we must make the conversion to physical imits of GeV. In j4l!] a 
continuum extrapolated value for the ratio v^/ S'^ was given for several gauge groups; in particular, ^/a^/ w 0.3351 
for SLr(2). Further, ^/a stands for the square root of the string tension. For this quantity, we used the input value of 
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FIG. 3: The optimal cr(fc^) in function of in units of = 1. 
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FIG. 4: The optimal k^Q (fc^) in function of k^ in units of GeV. The lattice (our analytical) results are indicated with triangles 
(dots). The error bars on the lattice data are roughly of the size of the triangles. 



^/a = 0.44 GeV as in (?]. Therefore, for SU{2), we find, 

2 

w 0.0194GeV2. (78) 

In FIG.m we have plotted the lattice as well our analytical result for the ghost dressing function, k^Q{k^), in units 
of GeV. We used the numerical data of fl^ l42ll , adapted to our needs. We observe that for sufficiently large k^, the 
lattice data and our analytical results converge, hi this case, the novel mass becomes zero as advocated earlier, 
meaning that we are back in the usual Gribov-Zwanziger scenario. For smaller k^, we found it more instructive to 
compare the lattice estimate of u(fc^) with our value as the errors on k^Q{k^) = becoming large when 

looking at a{k^) close to 1. FIG. |5] displays a{k^) in units of ^ = 1 up to fc^ = l x (^^^ = 0.0194 GeV^. We see 

that both results are in reasonable agreement, especially if we keep in mind that we have only calculated cr{k^) in a 
first order approximation. 
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FIG. 5: The optimal cr(A:^) in function of in units of = 1. The lattice (our analytical) results are indicated with triangles 

(dots). 



C. The gluon propagator 

We can apply an analogous procedure for the gluon propagator. This propagator at zero momentum, X'(O) is 



displayed in FlG.[6l We immediately see that there is an extremum at Ivfi = 0.33 y^^j / resulting in 

0.24 



V{0) 



(79) 



Domg the conversion to physical units again, we find for SU{2), 

12 

2?(0) ^ 

GeV^ 



(80) 



This value can be compared with the bounds derived from a partially numerical and partially analytical derivation 
Gl: 



1.2 12 

2 <^ < 2 • 



GeV 



GeV^ 



(81) 



We recall that our value X'(O) = 12/ GeV^ is a first order approximation and is only of qualitative nature. Neverthe- 
less, this value is still consistent with the boundaries of the lattice data set in fT3]. 



D. Violation of positivity 

We shall mvestigate if a gluon propagator of the tjrpe l (50|l displays a violation of positivity, another fact which is 
reported by the lattice data [2]. Following the analysis of [2J, the Euclidean gluon propagator can be expressed by a 
Kallen-Lehmann representation as 



+ O0 



dM- 



2 P(M^) 

^2 + M2 ' 



(82) 



whereby the spectral density p{M-^) should be positive to make possible the interpretation of the fields in term of 
stable particles. One can define the temporal correlator 0] 



C{t)= / dMp{M^)e 
Jo 



-Mt 



(83) 
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FIG. 7: C{t, IVfi) for a few values of t in function of M^, in units of fm. 

which is certainly positive for positive p{M^). However, if C{t) becomes negative for certain f, p{M^) cannot be 
positive for all M^, indicating that the gluon is not a stable physical excitation. Since C(f) can be rewritten as 

1 /'+°° 

C{t) = — e-^P'V{p^)dp . (84) 

2TZ J -co 

we must in fact calculate the ID Fourier transformation oi T>{p^). In FIG. [7] the Fourier transforms, C(t, M^) are 
shown for different t in units of fm. 

To determine M^, we again rely on the variational setup. We observe that for small t, there is no real extremum. 
However, for a certain t ~ 1, an extremum emerges at ^ 0. This extremum starts to grow for increasing t, but at 
the same time, the curve flattens out. Therefore, starting from t ^ 2.6, the extremum disappears again. Hence, we 
have taken'' = for t < 1, and = 0.18 for t > 2.6. The resulting temporal correlator is displayed in FIG.|8l 



3 jV4^ ~ 0.18 is the maximal extremal value, corresponding to f ~ 2.6. 
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We clearly observe a violation of positivity, which is in agreement with the lattice data (see FIG. 4 of Although 



C(t) (fm) 
0.5 r 



-0.1 
-0.2 



* • 4 • * 



t(fm) 



FIG. 8: C{t) in terms of t in units of fm in the refined GZ case, 
this agreement is only at a qualitative level, the shape of C{t) is very similar. 

It would be interesting to have a look at the temporal correlator in the pure Gribov-Zwanziger case {M^ = 0). 
Therefore, C{t) is displayed in FIG. |9] We can conclude that this plot is grosso modo the same as in the refined 
Gribov-Zwanziger setup. 

VII. CONCLUSION 



In this paper, we have extended the 4D analysis of the Gribov-Zwanziger scenario to 3D. The motivations for 
our work were to include an additional nonperturbative contribution to the Gribov-Zwanziger action, related to a 
dynamical mass generation effect, and the recent lattice data of il2ill4l1 . Two of the most striking features of these 
simulations were a finite nonzero value of the gluon propagator at zero momentum and the non-enhancement of 
the ghost, in both 3D and 4D cases. 



We have started from the following renormalizable action. 



Sgz + M2 / d^x 



a;«a;«)+2 



1) 2 
G7 



(85) 



which is an extension of the ordinary Gribov-Zwanziger action Sgz- Next to the local composite operator Ipcp — coco 
we have also added an extra vacuum term, necessary to ensure the restriction to the Gribov region O. In addition, 
the two parameters 7 and g are determined by the following gap equations. 



dY 



da{Q) 



aM2 



. 



M2=0 



(86) 



Before fixing M^, we have written down the explicit expression for the ghost and the gluon propagator. Firstly, 
the one loop ghost propagator Q (k^) is given by. 



fc2 1 - (7(fc2) ' 

with (7(fc2) given by expression ||63] |. At zero momentum, (7(0) is given by 



a{0) 



(87) 



(88) 



already after applying the gap equations. We see that for 7^ 0, cr{0) is smaller than 1, resulting in a non-enhanced 
ghost propagator. Secondly, the tree level gluon propagator T>{p^) yields. 



p4 _^ ^2p2 + ^4 ' 



(89) 
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FIG. 9: C{t) in terms of t in units of fm in the pure GZ case. 

with = Ig^Nj'^. Therefore, already at tree level, we find that 25(0) = which is nonzero for 7^ 0. We have 
also presented the one loop gluon propagator at zero momentum in equation | [60l l. 

In the last part of this paper, we have discussed a variational method in order to obtain an explicit value for in 
the Gribov-Zwanziger model, which has provided a value to the zero momentum ghost propagator (or equivalently 
(7(0)) and gluon propagator. With this method, we have found 

(7(0) = 0.94 , V(0) = . (90) 

Both values are in qualitative agreement with the lattice data. With this variational method, we also demonstrated 
the positivity violation of the gluon propagator, which is also confirmed by lattice data. 

Let us also spend a few words on the applicability of a weak coupling expansion. The reader will appreciate that 
we are no longer perturbing around a trivial vacuum, but are in fact considering a perturbative expansion around a 
nonperturbative vacuum characterized by a nonvanishing Gribov mass A, and additionally 7^ in our refined 
framework. These parameters ensure that the no-pole condition l|44l l is fulfilled, which puts a nonperturbative 
restriction on the theory. Obviously, we do not claim that we have all the relevant nonperturbative dynamics 
enclosed in this formalism, but at least the results we have do qualitatively match the available lattice data. The 
validity of a systematic perturbative expansion around this vacuum state is reflected in a coupling constant which 
should be sufficiently small, see e.g. | |5H] |. This reasoning also applies to the 4D case, see llStl . 

In summary, we have presented a framework, which consistently accoimts for the recent large volume lattice data 
in the infrared region in 3D and AD. The question is what happens in 2D, as the most recent data in 2D keeps 
predicting an enhanced ghost in combination with a vanishing zero momentum gluon propagator, contrasting with 
the higher-dimensional cases ||T3, fisl, lisll. Details of the Gribov-Zwanziger framework in 2D shall be presented 
elsewhere, as the situation is rather different there (43]. 

Previous studies on the infrared behavior of the gluon and ghost propagators within the Schwinger-Dyson 
formalism in 4D which where in compliance with the lattice data can be foimd in llsl l46l |4^ . In particular, we 

observe that in [4^ a gluon propagator fit was given as f (p^) = 4 ^ ^"1° — 4, while the ghost propagator almost 

behaved like It is interesting to notice that this kind of gluon propagator is of the same type as the one found 

here, in a completely different way. Finally, let us also mention that similar research in the maximal Abelian gauge 
| |48|| also gave results in qualitative agreement with the available lattice data in that gauge, see 1 , 49,] . 

VIII. ACKNOWLEDGMENTS. 

We are grateful to A. Cucchieri and T. Mendes for providing us with lattice data and useful comments. D. Dudal 
and N. Vandersickel are supported by the Research Foundation - Flanders (FWO). The Conselho Nacional de Desen- 
volvimento Cientifico e Tecnologico (CNPq-Brazil), the Faperj, Fundagao de Amparo a Pesquisa do Estado do Rio de 
Janeiro, the SR2-UERJ and the Coordena(;ao de Aperfei(;oamento de Pessoal de Nivel Superior (CAPES) are grate- 




t(fm) 



22 



fully acknowledged for financial support. This work is supported in part by funds provided by the US Department 
of Energy (DOE) under cooperative research agreement DEFG02-05ER41360. 



APPENDIX A: CONSTRUCTION OF THE EFFECTIVE ACTION T{a) FROM THE GENERATING FUNCTIONAL W(/) 

1. General framework 

We define the classical field cr{x) conjugate to the source J{x) as follows^: 

with W(/) the generating fimctional defined in l fT6l l. The effective action T{a) is then obtained in the usual way by a 
Legendre transformation 

Tia) = W{J) - J d'^xJix)o-ix) , (A2) 
where }{x) is understood to be a functional of a{x). It is easily derived that 



S(7{x) 



}{x) . (A3) 



The original theory is recovered when the source / attains again a zero value. A source is nothing more than a tool 
to probe the theory. In our case, we are incorporating the effects induced by the operator fcp — coco. At the end, 
a source should always become zero. The equation iA3i for example corresponds to the quantum version of the 
classical equation of motion when / = 0. 

In what follows, we shall limit ourselves to constant / and cr, as we are mainly interested in the (space time 
independent) vacuum expectation value of the operator coupled to the source /. The remaining task is to obtain 
the functional form of F in terms of a. Usually, this is done using the background field formalism of Jackiw 
150], when an effective potential T{(pc) associated with an elementary field ^(x) is determined, whereby it is 
understood that (p{x) = (pc + (p{x). The quantity (p{x) represents the quantum fluctuations around the vacuum 
expectation value of(p{x), i.e. {(p{x)) = (pc, {(p{x)) = 0. Unfortunately, when a composite operator is considered, this 
procedure is less useful, as there is no elementary field associated to the operator to begin with. Sometimes, tech- 
niques are at one's disposal to introduce such a field, e.g. by employing a Hubbard-Stratonovich transformation [2^. 

Let us proceed thus by explicitly performing the Legendre transformation, along the lines of ll5lll52ll . The generat- 
ing fimctional W(/) will be obtained as a series in the coupling constant, 

W(/) = Wo (/)+ ^2^1 (/) + ... . (A4) 

As a consequence, we may write 

^(/) = ^^0(7) +gVi(/) + . . . , with anU) = ^ . (A5) 
This last series can be inverted to give / as a function of cr, 

/(^) = h{<T)+g^h{<T) + ... . (A6) 

As a trivial consequence, 

C7 = (7(/(C7)) = a{Jo)+g^hcr^i}o)+a^{h))+... . (A7) 



* To avoid any confusion with the reader, this classical field a{x) has obviously nothing to do with the one loop correction to the ghost form 
factor i7(fc^) first defined in equation iiol . 
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Since the field a is supposed to be independent of the coupling constant at the level of the inversion^, we find an 
iterative inversion procedure, 

= o-qUo) =^ Jo = /o(t^o) / 
^i(/o) 



= /ic7-^(/o) + ^i(/o) h 



(A8) 



In this fashion, every term in the series for / can be expressed in terms of Jo, which itself is a function of ctq. 
Summarizing, everything can be written in terms of ctq. Doing so, we can calculate the r.h.s. of l |A2| l up to the desired 
order in by simply substituting the corresponding expressions for /,■ in the l.h.s.. 

Once the inversion is performed, we must fix cq, which corresponds to the condensate of the operator coupled to 
/, by demanding that 

=0^/0(^0)+ M^o)+...=o. (A9) 

dCTQ 

2. Application to the LCO Tpq) — cocv 

We shall now employ the previous method on the 3D Gribov-Zwanziger action, in the presence of the LCO 
'^cp — coco. The action we use is given by (|26l l. Using (|A5l l and llASt , yields 

[3 ^ 3 /^/o + v//o-4A4^ }, \ 3//o-V/o-4A^W. h 




(AlO) 

Defining a = j^^ctq, we found 6 different branches for the inverse function }o{c)- However, only one of these was 
real valued and gave rise to a solution of the gap equation, so we focus our attention on this particular solution. 



/0('^) = ^2 + 



3 ;v 



8192^8 + 442368A2c^ + 6469632A*?* + 20901888Af'52 + 45349632A8 + 186624 v^A5(9A2 + 4P)\/243A2 + 8?^ 



27^2 432^22^2 

256c^** + 9216A2cf*' + 51840A*^* 



21622/3^2 3/2^4 |^8192c^8 +442368A25^6 +6469632A4?4 +20901888A%2 _|_45349632A8 + 186624 v^3A5{9A2 + 4cf2) V243A2 + 8?2j 

We have displayed a plot of /o/A^ in terms of (7 /A in FIGHO] Solving the gap equation Joicr) = numerically leads 
to 

- 1 n 

ao w 1.06— A w 0.056(n2 - 1)A . (All) 

on 

This numerical solution is consistent with the analytically derived solution l|22] |. Let us now determine the Gribov 
mass A. Therefore, we first compute the effective action using its definition l|A2|l . We did not write down the eventual 
result in terms of A and a, as the expression is quite lengthy. We subsequently solved the gap equation (horizon 
condition) = (see FIG ITT) numerically in the case N = 3, and we found 

A w 0.113/ , (A12) 



which is also consistent with the analytical result ||24j. The corresponding vacuum energy is given by 

Evac « 0.000214/ , (A13) 



^ Only after that the gap equation, ^ = 0, is imposed, a will collect a g^-dependence. 



24 



4 - 



3 - 




1 2 3 4 5 a 

FIG. 10: A plot of ^ in terms of |. 
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FIG. 11: The horizon function |j for N = 3 in units = \. 



again equivalent with ||25)|. 

The gap equation derived from the effective action r(£7-o. A) should of course give back this perturbative solution, 
next to a potential nonperturbative solution, which can never be discovered by simply using JlSt . The whole point 
of the inversion (Legendre transformation) is just that there might be multiple solutions to the equation / = 0, next 
to the perturbative one. hi our case, the situation is interesting because there is already a nontrivial scale in the 
game, namely the Gribov mass A. This allows us to obtain a nonvanishing value for the condensate already at the 
perturbative level. 

APPENDIX B: THE RENORMALIZABILITY OF THE EXTENDED GRIBOV-ZWANZIGER ACTION: ALGEBRAIC 

ANALYSIS 

1. The Ward identities 

The action l l27l l exhibits several Ward identities which we have enlisted below. 
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The global U{f) invariance: 

Uiji: = 0, (Bi) 



The Slavnov-Taylor identity: 



5(E) =0, (B2) 



5(E) = d''xl—— + —— + b"— + w"— + Lo"— + M'!! r + N', 



SKf, SAf, SL" Sc" Sc" 6TD1 ' 6(p1 '' SUf,' ^ 5Vf 



• The Landau gauge condition and the antighost equation: 



The ghost Ward identity: 

Q"^ = A«i, (B5) 



Notice that A"j is a classical breaking, since it is linear in the quantum fields. 
The linearly broken local constraints: 







SMfl 

r 






ST. 












SI. 




sui; 


SI. 




SI. 



gf^AlV^Z + Jcpl, (B6) 



S¥ 

SI. 



^-=gf''^'A';u;' + Ja>^ , (B7) 



The exact 72.,^ symmetry: 

Uijl. = 0, (BIO) 



^ m ' sn:j ' '^M- 



We shall now rely on the algebraic renormalization procedure 113111 , according to which the most general allowed 
invariant counterterm E*^ is an integrated local polynomial in the fields and sources with dimension bounded by 
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three, with vanishing ghost number and Q yr-charge, and subject to the following constraints: 



SI' ^ SI' rabc-M^^' 

SI' ^ SI' 



suj"' ' sm 



Uiil' 


= 0, 




n 


SZ SI 

''SKf, 


= 0, 


SI' 
'Sb" 


= 0, 




= 0, 


xyc 

„ fabcj jbi 




rabcybi ^'^'^ 

^ SKI, 




SYf ^ SI' 

: + du : 


= 0, 


slo"' ' sm 


SI ^ sz 

: + 3,, r 


= 0, 


S-^"' ' SMf/ 


TZijl' 


= 0. 



(Bll) 



The operator Bj; is the nilpotent linearized Slavnov-Taylor operator, 

„ [ ,^ ( SI S SI S SI S SH S ,^S _^ S „ S s s 



SK"^, SAf, SAf, SK^ SL'' Sc" Sc" SL'' Sc" ScJf ' S<p1 ^ SUf ^ SVj'/ I ' 
ByJB^ = 0. (B12) 

One can show that HF does not depend on the Lagrange multiplier W , and that the antighost c" and the f- valued 
fields (fl, (jjI, oJ" can only enter through the combinations 

rnbci^jbi ,„ci „rabc-trbi—ci 



n;' = n;; + a,,a;«' , m;; = v;' + a,,^"' . (bb) 

Imposing the previous constraints arising from the Ward identities, the most general counterterm can be brought in 
the following compact form 

1' = aoSYM + a,j d'x a; ^ + K^d^c" + t/^M- - l1,TN«'^ + j d^xpg^} , (B14) 

with flQ/ 1^1 arid fl2 three arbitrary parameters. It is an impressive feature of the action Z, that only 3 independent 
parameters can enter the coimterterm, despite the presence of many fields and sources. 



2. Stability of the action and the renormalization constants. 

As the most general counterterm 1' compatible with the Ward identities has now been constructed, it still remains 
to be checked whether the starting action E is stable, i.e. that 1' can be reabsorbed into E through a renormalization 
of the parameters, fields and sources appearing in E. 

According to expression (|B14b , 1' contains three parameters aq/^i arid a2, which correspond to a multiplicative 
renormalization of the gauge coupling constant g, the parameter p, the fields (p = (A",, c", c", b", cp^, col, '^1' '^1^ 
the sources O = (X"?', L" , Mf, Nf, V^f', !J«', /), according to 



l{g,p,,p,<^) + r}i:' = I{go,po,cl>o,^o) + O{r]^) 



(B15) 
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This is possible, provided that we define the bare quantities in terms of the renormalized quantities as 

^0 = Zgg , po = Zpp , 

(l>o = zy^cp , Oo = Zo* , (B16) 

with 

Zs = 1+Vy' 2y' = l+'?(«i-y) ' Zp = l+ri{a2-a^-ao) . (B17) 

These are the only independent renormalization constants. For the rest of the fields, we have 

Zc = = Z(p = Z-^ = Zi_Q = = Zg Z^ , (B18) 

while for the renormalization of the sources (^M"/, V^f , W^}, 

Zm = Zjv = Z,. = Zu = Z^-i/22^1/4 ^ 2j = Z^^zy2 . (b19) 

We see thus that the LCO Ipcf — To to does not renormalize independently, as it is evident from l|B19t . The only new 
parameter entering the game corresponds to the renormalization of the vacuum functional, expressed by p and its 
renormalization factor Zp. As discussed in the main body of this paper, this parameter turns out to be equal to zero 
anyhow. 



APPENDIX C: PROPAGATORS IN LATTICE AND CONTINUUM FORMULATION 

For the benefit of the reader, in this Appendix we discuss how our results compare with the corresponding lattice 
data in the case of the gluon propagator. The quantity which is evaluated in both cases is the gluon propagator, 

namely the connected gluon two-point function where the gauge field configurations A", are re- 

stricted to the Gribov region Q. We shall also show that the gluon and the ghost propagators are color diagonal. 



1. Continuum fonnulation 



In the continuum formulation, the gluon propagator is given by the connected gluon two-point fimction, and 
expressible by means of 



A«(x)A^(y) 



(CI) 



/=o 



with Z'^(/) the generating functional of the connected gluon n-point functions, which in our case will read 

-z^{/) ^ = |[dY]e'(''°'+-^'''K^/"^?'+^^''<+^l!;'0), (C2) 



where Stot is the improved Gribov-Zwanziger action | [46l l. This amounts to considering the Landau gauge fixing, 
such that the relevant gluon configurations belong to the Gribov region, i.e. these are (local) minima of / d?xA^ 
along the gauge orbit. 

As proven in lITsIl , the gluon propagator (|Cll l is transverse. In a condensed notation, one shall find 



( ^ \ 

} }f }f + higher order terms in /, /<p, /-^ 



(C3) 



where A4 is the matrix propagator, as written down in | [55t up to first order. The upper left corner of this matrix A4 
corresponds precisely to the gluon propagator, as it is apparent by taking the second derivative with respect to the 
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source /J^. 

If one is interested in the IPI two-point functions, one should look at the corresponding generator, which is 
the effective action r[A^,, cp,f]. As is well known, the corresponding IPl two-point function will be the inverse of 
the connected two-point function (or propagator). Said otherwise, the corresponding matrices will be each others 
inverse. This is also explained in the main body of the text, with the IPl two-point function matrix written down in 
(|53] l, again up to first order. 

As we have already stressed earlier in the paper, the extra fields (^"^ , cp^f ,cJlf , lo"I^^ are introduced in order 

to obtain a local manageable field theory, which is capable of restricting the gauge field configurations to the Gribov 
region, which is a rather nontrivial operation in the continuum. In principle, one could opt to work in an effective 
field theory fashion by again integrating out the extra fields. Clearly, this will give rise to a very complicated 
(nonlocal) action, written solely in terms of the original Yang-Mills fields. In this formulation, the gluon propagator 
is directly related to the inverse of the IPl two-point function, due to the absence of mixing. Anyhow, the result 
for the propagator itself will be the same as the one already obtained before in the preferable local and manifestly 
renormalizable formulation with the extra fields, when looking at the same order in g-^. This can be easily checked 
at tree level: integrating out the auxiliary fields in l l46t leads to the following quadratic (nonlocal) effective action. 



■Squad / d " 



(C4) 



where the limit a — > is understood in order to recover the Landau gauge, and where we skipped the irrelevant 
constant terms. The tree level gluon propagator in momentum space is in this case 1 / Q2, with Q2 the quadratic form 
appearing in ||C41 , when expressed in momentum space. Clearly, this leads back to the lowest order approximation 
in the upper left corner of (|55] l. 

Concerning the ghost propagator, a similar formalism applies. 

Let us discuss the issue of color diagonality. The global color symmetry guarantees us that the gluon and the 
ghost propagators are color diagonal. This property is encoded in the global SU{N) Ward identity, which reads at 
the classical level 

with Z the classical action and cp all the other fields. In particular, 

^U^; = f'^i, , (C6) 

and similar relations for the other fields (p. Therefore, we find 
This identity can be upgraded to the quantum level^, 

with r the generator of IPl correlators. Performing the Legendre transformation leads to the analogous Ward identity 
for the generator Z'^ of connected correlators, 

) = 0. (C9) 



6 We refer to [31] for the explicit proof. 
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We shall concentrate on the gluon sector, so we have already set all other sources equal to zero. Taking a derivative 
w.r.t. leads to 



r 



ibc 



. 



(CIO) 



Next, taking a derivative w.r.t. /^(z) and setting / = at the end gives the following relationship. 



■die 



S^Z" 



5'^Z' 



6rAz)6]i{y) 



0, 



(Cll) 



or equivalently 

f'^{Al{y)A{{z))+f''^{Ai{y)A\{z)) = 0. (C12) 
This relation expresses nothing else than that the gluon propagator is an S{J(N) invariant rank two tensor. Therefore, 

(A^,(y)Ai(z)) a 5^', (C13) 
since 5^^ is the unique invariant rank two tensor. 

Obviously, all available explicit loop results obtained with the (modified) Gribov-Zwanziger action are com- 
patible with the general proof. Notice also that the proof is the same as the one we would use in the case of normal 
SU{N) gauge theories. 



An analogous result can be derived for the ghost propagator. 



2. Lattice formulation 



hi a lattice formulation, one also calculates the connected two-point function, by taking the Monte Carlo average 
of the discrete version of the operator (^A1,{x)Ay{y)'^. The statistical weight for this simulation is given by the 

exponential of the discretized version of the Yang-Mills gauge action, e.g. the Wilson action. The Landau gauge 
fixing is numerically implemented by minimizing a suitable functional along the gauge orbits, which corresponds 
to minimizing J d^xA^ in the continuum. As we have already explained in the introduction, this amounts to nu- 
merically selecting a gauge configuration within the Gribov region, equivalent with what we did in the continuum. 
We refer the interested reader to Section 2 of [53] for the explicit expressions of the discrete action, gauge fields and 
minimizing functional. In particular, we refer to Subsection 2.4 in which the continuum and lattice versions of the 
gluon propagator are written down. The lattice gluon propagator also turns out to be transverse. Moreover, both 
the gluon and ghost propagator are found to be color diagonal. 

We emphasize here that lattice simulations thus never directly calculate any IPI two-point function, but, we 
repeat, do also calculate the (connected) two-point correlator, i.e. the propagator itself. 
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